We propose highly accurate finite-difference schemes for simulating wavepropagation problems described by linear second-order hyperbolic equations. Theschemes are based on the summation by parts (SBP) approach modified forapplications with violation of input data smoothness. In particular, we deriveand implement stable schemes for solving elastodynamic anisotropic problemsdescribed by the Navier wave equation in complex geometry. To enhance potentialof the method, we use a general type of coordinate transformation andmultiblock grids. We also show that the conventional spectral element method(SEM) can be treated as the multiblock finite-difference method whose blocksare the SEM cells with SBP operators on GLL grid.
展开▼